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We consider quintessence scalar field cosmology in which the Lagrangian of the scalar field is 
modified by the Generalized Uncertainty Principle. We show that the perturbation terms which 
arise from the deformed algebra are equivalent with the existence of a second scalar field, where the 
two fields interact in the kinetic part. Moreover, we consider a spatially flat Friedmann-Lemaitre- 
Robertson-Walker spacetime (FLRW), and we derive the gravitational field equations. We show 
that the modified equation of state parameter wgup can cross the phantom divide line; that is 
wgup < —1. Furthermore, we derive the field equations in the dimensionless parameters, the 
dynamical system which arises is a singular perturbation system in which we study the existence 
of the fixed points in the slow manifold. Finally, we perform numerical simulations for some well 
known models and we show that for these models with the specific initial conditions, the parameter 
Wgup crosses the phantom barrier. 
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1. INTRODUCTION 

Various approaches of quantum gravity, such as. String Theory (ST), Doubly Special Relativity (DSR), Black 
Hole Physics (BHP) predict that there should exist a minimum length scale of the order of the Planck length (Ipl), 
or, in other words, the existence of a maximum energy scale in nature. This prediction motivates to modify the 
Heisenberg’s Uncertainty Principle to some Generalized Uncertainty principle (GUP) in the context of the quantum 
theory of gravity [l|. The GUP model that has been proposed from String Theory and Black Hole Physics contains 
the existence of minimum measurable length, whereas DSR provides a GUP model which incorporates the existence of 
a maximum observable momentum. This maximum momentum scale has been obtained in the context of DSR theory 
for the consideration of the frame independent invariant upper bound of energy scale into the theory . Interestingly, 
the above two scenarios, the minimum length scale and the maximum momentum scale can be combined into a single 
deformed uncertainty principle (DUP), which has the dual nature Q. However, the modified Heisenberg’s algebra that 
can be generated from this DUP relation is not relativistically covariant. Hence, it is very difficult to apply the DUP 
in any relativistic scenario. So, here we mainly concentrate on the relativistically covariant modified Heisenberg’s 
algebra, which is constructed from the GUP model that incorporates the minimum length scale only. Indeed, later 
on, we will see that this modified algebra deforms the coordinate representation of the momentum operator, whereas 
the position operator remains undeformed. 

On the other hand, standard General Relativity (GR) is not enough to describe the phenomena that recently 
appeared in astrophysics and cosmology, such as, the late-time acceleration of the universe [MU. An easy way to 
explain the accelerated expansion of the universe is to consider an additional fluid which has a negative equation of 
state parameter (dark energy); this new fluid counteracts the gravitational force, and leads to the observed accelerated 

S ansion. The simplest dark energy candidate is the cosmological constant A leading to the ACDM cosmology (l^ 
However, this model has two serious problems: the fine tuning problem [13)[13, and the coincidence problem [la| . 
As a consequence, different dark energy models gave rise into existence, and try to explain the recent accelerating phase 
of the universe, such as, time-varying A (t) cosmologies, quintessence, k-essence, phantom and others (for instance see 

Pil)- 
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Now, it is known from cosmic microwave background observations, that our universe had gone through the in¬ 
flationary phase at its early stage, and then entered into the radiation, and matter dominated eras, in which the 
expansion became decelerating. So, at the early stage of high energy regime, it was accelerating, and, in this present 
low energy regime, it is again accelerating. Quantum gravity effects played an important role at the early universe, 
see for instance [241 - 1^ . For instance, the mystery behind the origin of magnetic fields with /iG (microgauss) strength 
(having approximately the same energy density as the cosmic microwave background radiation) observed in 
the intergalactic scales, has some explanation in the context of GUP [s^ . Furthermore, the existence of a minimum 
length, specifically GUP, prevent the complete evaporation of black holes [sj . In particular it has been shown that 
a black hole with temperature greater than the ambient temperature should radiate photons, as well as, ordinary 
particles until it reaches the Planck size, whereas in the Planck scale, it stops radiating, and, its entropy reaches zero 
[S^ (for a review see [s^). 

Thus, it will be very interesting to investigate whether the GUP corrected modified theory can explain the present 
day accelerating phase, or, not. Because, these theoretical approaches predict that the minimum length effects are 
quite universal at any stage. Hence, they must have to describe the late time accelerating phase in order to keep the 
consistency in the theory. Since, it is known that the dark energy, is the proposed candidate to explain the present 
day acceleration; so, we can study the consequences of the GUP modification on the dark energy models. 

We have already mentioned that the modified Heisenberg’s algebra corresponding to the GUP model that in¬ 
corporates the existence of minimum length, modifies the coordinate representation of momentum operator. As a 
consequence, the dispersion relation of momentum become modified and finally the Klein-Gordon equation deforms 
accordingly. From this modified Klein-Gordon equation one can easily expect the deformation of a Lagrangian density 
in field theory. As, in the context of dark energy we have a field theoretic model, namely, the quintessence scalar field 
models ^ [O [l3. ISGI - I^ . so, we can modify it in light of GUP. Therefore, in order to study the late-time acceleration 
of the universe, in this work, we are interested on the scalar field models (quintessence) in which the Lagrangian has 
been modified by the GUP. Note that, for this Lagrangian the gravitation equations follows from the Einstein-Hilbert 
action with the usual variational principle. 

From the analysis of the recent data @-@1 it appears that the EoS parameter for the dark energy could have values 
wde ^ ~1- The limit wde = is called phantom divide line or phantom barrier. As, the usual quintessence scalar 
field model can not explain this phase, so, our another motivation is to modify the scalar field model in such a way 
that it can explain the recent phase. Indeed, here we will show the modified model will drag a correction term into 
the effective equation of state parameter wgup ( 2 ), and it concludes that, wgup {z) can cross the above phantom 
barrier. The plan of the paper is as follows. 

In section[2l we give the basic definitions and properties of GUP and we derive the modified Klein-Gordon equation. 
In section [31 we consider a quintessence scalar field cosmology follows from GUP. We use Lagrange multiplier and 
we write the gravitational action as that of two scalar field cosmology with a mixed kinetic term. Furthermore, we 
consider a spatially flat Friedmann-Lemaitre-Robertson-Walker (FLRW) spacetime and write the field equations and 
the modified EoS. In section |4l we write the field equations in the dimensionless variables, the dynamical system 
which arises is a slow-fast dynamical system, where we study the fixed points of the system in the slow manifold. 
Furthermore, we perform numerical simulations for two well known scalar field potentials. Finally, in section [51 we 
draw our conclusions. 


2. GUP-BASED MODIFIED MODEL 

The GUP consistent with the existence of the minimum measurable length has the structural form as follows 

AX^APj > + I3P^) + 2pP,Pj]. ( 1 ) 

The deformed Heisenberg algebra that can be found from the above GUP model H]) is given by , 

[A„ Pj] = ih[S,,il + /3P2) + 2/3P,P,]. (2) 

Here /3 is a parameter of deformation defined by ^ = Po/MpiC^ = where Mpi is the Planck mass, ipi 

10“^^ m) is the Planck length, Mpic^ 1.2 x 10^^ GeV) is the Planck energy. 


^ A quintessence scalar field has a time varying equation of state (EoS) parameter Wfj) {z) which is bounded as follows {z)\ < 1, and 
could describe the acceleration of the universe when (z) < — 
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From the commutation relation ([2]), the coordinates representation of the momentum operator can be modified to be 
Pi = Pi(l + /3p^)) by keeping = Xi undeformed. Here {x,p) is the canonical representation satisfying [xi,pj] = ihSij. 
In the relativistic four vector form, the above commutation relation ([2]) can be written as 1^ 

[X^,P,] = -ih[il - - 2f3PM (3) 

where the flat metric is given by = (1,-1,—1,-1). The corresponding deformed operators in this case are 

PfJ. = PaP'i)) , = x^, (4) 


where p^ = and [Xfj_,p„] = 

For a spin-0 particle with rest mass m, the Klein-Gordon equation is p^^PfiPv^ — (mc)^ 'k = 0. By substituting 
Pfj, from (ID), we have the modified Klein Gordon equation has the following form 

+ 2phA (A^-) -h % -k O (P^) = 0, (5) 

where A = [d ^''is the Laplace operator, whereas for the metric —>■ Vfj.u, A = □ = 9^9^. In 

contrary with the classical Klein-Gordon equation in equation ([S|), there is the term 2/3/iA (A'k), which is the quantum 
correction term. Furthermore, we remark that equation ((S]) is a fourth-order singular perturbed equation. 

Equation JS]) could arise from a variational principle, that is, there exists a function L ...), such that 

E (L) = 0, where E is the Euler Lagrange operator. The Lagrangian of the fourth-order Klein-Gordon equation has 
the following expression (i^ . 


where the new operator is 22^ = -I- (A); is the covariant derivative, i.e. and e = pfP = 

Pdlil2. 


3. GRAVITATIONAL ACTION AND FIELD EQUATIONS 

Let a four dimensional Riemannian manifold with metric and curvature i?. We consider the following gravita¬ 
tional action 


S = Sgr + Soup + 5 ' m , ( 7 ) 

where Sqr = / dPx^—gR is the Einstein-Hilbert action, Sm describes the matter components, and Soup is the 
action of the scalar field modified by the GUP, that is, 

Soup = -J^ (^^g^’^V^pVyP - . (8) 

Therefore, the gravitational action is^: 

S = Jy/^d*x(^R- ]^g^''V^4>Vy(t> + V{4>)^ + Sm , 

or, equivalently, 

s = Sgr + s^ + SmP J {-eg>^''(j).(^fj,g°‘''(j).\ab\^)), (9) 


^ In contrary to the Saez-Ballester scalar field model in which the operator T)^, can be seen as / (0) where is the covariant 
derivative, the new operator includes second-order derivatives. For the cosmological viability of the Saez-Ballester model see [43 . 
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where is the action of the “classical” scalar field [1^. Note that, the last term of the above action can be 
written equivalently as follows 


'5'c = y ■ (10) 

The action integral ([9]) provide us with a higher-order theory of gravity in which the higher-order terms follow from 
the term Sc- 

In order to simplify the calculations, we introduce a Lagrange multiplier^, and, consider the new variable ip defined 
hy Ip = g^'''(p-^^v\ therefore, incorporating the Lagrange multiplier A into (fTOll . and by solving the Lagrange’s equation 
^ = 0, we have A = —2eip. Hence, applying integration by parts, equation (nni becomes 

'S'c = y d'^x (-2eg“V,aV'.fc - £'0^) ■ (H) 

Therefore, the gravitational action (|9l) becomes 

S = SgR + + Sij; S Sm, (12) 


in which. 


>5'v> = y d'^x [-2eg°‘^(p,aip,b - £0^] 

The last action describes a field which is coupled with the field (p, where the extra term follows from the perturbation 
effects of GUP. 

Finally, from action (USD, we have that the gravitation equations are 

= R^v - -Rg^y = 8TrGTf_ii, (0,0) -I- S-nGT^^ , (13) 

where {(p,'^) is the energy-momentum tensor for the two fields'* {0,0}, and, is the energy- momentum of the 
matter component. We remark that the kinetic metric which is defined by the two fields {0,0} can not be written in a 
diagonal form, since the “coordinate” transformation {0,0} —> {0,0} which diagonalizes the kinetic metric depends 
on £“*, i.e. it is singular transformation. Which means that the action (1121) can not be written in that of the Quintom 
model 1^. Furthermore, the second fluid 0, is not a real fluid but it follows from the higher-order derivatives of 
the field equations. 

In particular, the reduction of the order of the field equations from a fourth-order theory to a second-order, by using 
a Lagrange multiplier, it means that in the same time the number of the dependent variables is increased. An other 
analogue is the / (i?)-gravity in the metric formalism [d^ . The last gravitational theory for nonlinear / (i?) function 
is reduced with the use of the Lagrange multiplier in a Brans-Dicke theory with vanished Brans-Dicke parameter , 
the so called O’Hanlon massive dilaton theory [5^. For other applications of Lagrange multipliers in different theories 
of gravity, see [sil - ls^ , and for a discussion on the degrees of freedom in alternative theories of gravity, see . 

In the following section we continue with the derivation of the field equations in a spatially flat FLRW spacetime. 


3.1. Field equations in a FLRW background 

Let us consider a spatially flat FLRW spacetime with line element 

ds^ = dt^ — (t) [dx"^ + dy^ + dz'^) , (14) 

and, with gravitational action m, where the matter tensor describes a dust fluid with energy density pm- 
Moreover, from the conservation law = 0 (Bianchi identity), and, for the space (fCT) . we have that 


3 

4 


Lagrange multipliers are useful to reduce to order of the differential equations. However at the same time, the dimension of the dependent 
variables is increased. For applications of the Lagrange multipliers in high order theories of gravity, see for instance 511453 
For exact solutions of two scalar field models, see [s^. and for a quintom scenario with mixed kinetic term, see [Hs 
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where pmo is the present value of the energy density. From the cosmological principle, we have that the fields {(f), il)'\ 
inherit the symmetries of the metric (HI, hence, </) = </) (t), 'ip = -tp {t). 

Therefore, for the space da, the gravitation equations (1131) follow from the following Lagrangian® 

L = 3ad^ — ~ ‘2:ea^^’ip + a^V {(p) — ea^'ijP' + pmo- (15) 

Lagrangian da can be seen as the Lagrangian of a particle moving in a Riemannian manifold with metric , and 

effective potential 14// = — [a^V (</>) — ea^tp^ + pmo), that is, 

L{y\y^) = - Veff (/) , (16) 

where, y* = (a, (p, ip), and yi/ is given by the following line element 

ds^ = Qada? — a^dcp'^ — Aea^dcpdip. (17) 

For the Euler-Lagrange equations of (I16|) we have that these are given by the following expression 

l^kf + l^rr^Jky^y’^ + (K//)., = o, (18) 


where F^,^ are the Christoffel symbols of jij, and the nonzero Christoffel symbols are 


and 


^ (1 F 

pa _ pa _ _ pa _ ^ 

^ ~ 2a ’ ^ ~ 4 ’ M ~ 


F"^ = — = — 

2a ’ 2a' 


(19) 


( 20 ) 


Therefore from da-da, we have that the second Friedmann’s equation which is the Euler-Lagragne equation of 
da with respect to the scale factor, a, is given as follows 


d + —+ ae^tp - -a [V {(p) - eip'^) = 0, 


( 21 ) 


and the Klein-Gordon equations for the “two field” {0,'*/'}, are 


3 • 

(j) H— a(j} — ^ = 0, (22) 

a 

which follow from the conservation law (</>, ip) = 0. Eurthermore, the first Friedmann’s equation can be seen as 

the Flamiltonian function which follows from Lagrangian da, and it is®, 

PmO = 3ad^ — 2 ®^^^ ~ 2ea^pip — a^V (cp) + aPeip"^. (24) 

Here, we would like to note that, the set of the field equations da-da form a singular perturbation dynamical 
system. Moreover, it is straightforward to prove that, equation (l22ll is the constrained equation ip = g^''(p-,fj.ij, and, 
if we replace ip from (1221) in (I23p . then we will have a fourth-order equation, which is the modified Klein-Gordon 
equation for the field, (p, in GUP. Furthermore, in the limit in which e = 0, that is, the deformation parameter /3 is 
vanished and we have the classical uncertainty principle, Lagrangian m is that of quintessence scalar field with dust 
fluid, and in the same time equations dUD-dMl) reduced to the field equations for the last model. 


^ In the following, we consider SttG = 1. 

® Alternatively equation |[24t can be derived from the Euler-Lagrange equation with respect to a new variable, which arises from the 
lapsed time dt = Ndr, for instance see [hsil . 
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An alternative way to write the field equations is, 

3 H^ = PGUP + Pm, (25) 

2H + 3H^ = -PGUP, (26) 


where pgup = P4, + PG, PGUP = P4>+ PG- Here, p^, p^ are the density and pressure of the usual scalar field given by 


p<i> = + y (^) 

, = -y (<(’), 

(27) 

and pc, pc are the perturbation terms 



pc = £ (^pTp - 'ip'^'^ 

, PC = £ (^p'tp + ■ 

(28) 


However, the two fluids {p^, pc} are interacting; furthermore, the field pc, i-e. equation (E51) is not a real field, but 
it was introduced by the quantum corrections of GUP. 

Finally, the effective EoS parameter for the scalar field which follows from GUP is 

+ e (^ 2 ^ 4 } + 

+ e (^2^ip — 

Using the Taylor series expansion of equation (1^ at e —>■ 0, we have wgup = w' , where = ^, is the EoS 
parameter for the “classical” scalar field and 

w' = — (pG - w^pc) + O (e^) . (30) 

P(j> 

In the limit where (jP V , which means that ~ —1, then w' = £y(j)'ip; therefore wgup becomes 

4e • • 

WGUP — —1 + -y4>4’ ^ when ipij} < 0. (31) 



WGUP = 






Hence, the EoS parameter wgup can cross the phantom barrier provided when (ptjj < 0, and, 




>> 1. Here, we 


would like to note that, since the field equations form a singular perturbation system (slow-fast system), it is possible 


that 


ip 


eip 


could be big enough in 


>> 1, which means that either when the parameter e is very small, the product 

order the perturbation effects to be measurable. In order to show this, in the following section we study the existence 
of fixed points for the singular perturbation system in the slow manifold. Einally, we perform numerical simulations 
for some well known scalar field models, and, we show that, wgup can cross the phantom divide line. 


4. COSMOLOGICAL EVOLUTION 


In this section, we consider the field equations (HI where now the energy momentum tensor T^i, is + 

where tIpP\ , are the energy-momentum tensors for a dust fluid {wm = 0) (dark matter), and a radiation 
fluid {wr = ^) respectively. 

The field equations (1131) for the spatially flat FLRW spacetime are 


3H^ — p^ + Pm + Pr + £ {^4’^ ~ F (,4’)^ , 
2H + =-p^ - Pr - £ (2pi; + F , 

where F {tp) = ip'^. 

We define the dimensionless variables 111,16^ 



VvW 

yj Pr 

Vqh ’ 

X2 — r- > 

V3H 

'^'6 ! — ? 

^/3i^ 

2\/6 ip 

Vf 


3 H 

’ V3H 

’ 3i72 


(32) 

(33) 


(34) 

(35) 
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In the new variables, the first Friedmann’s equation (I32[) can be written as follows 

— 1 ^GUP ^ 


(36) 


where now = ( 0 : 3 ) , and ^gup = (a^i) +( 2 ^ 2 ) +e — ( 2 / 2 ) 

Furthermore, the second Friedmann’s equation (155)) takes the form = |S (x, y,e), where 


2 / n 2 , / ^2 1 / x2 


S (x, y, e) = -1 - |^(a;i) - {X 2 ) J - e yxiyi + ( 7 / 2 ) J - ^ (^a) 

hence, the total EoS parameter is Weff = — 1 — S (x, y, e). 

The EoS parameter (1551) for the scalar field in the new variables becomes 


WGUP = 


(xi)^ - (X2)^ 

+ e 

^xiyi -f (y2)^) 

(Xi)^ -f {X2 f 

-f £ 

xiyi - ( 222 )^) 


where in the limit —>■ — 1, becomes 


WGUP = —1 + 2e--!^ + O (e^) . 
{X 2 ) 


(37) 


(38) 


The derivation of the variables x, y with respect to the lapse time N = In a gives the following dynamical system 

(39) 


^ = f (x, y. A, Ai, e) > ^ = *1 (x, y, A, e), 


^ = 52 (x, y. A,/i, e) , — =/i 2 (x, y. A, e), 

dyi r \ \ 

^d/V ^ (^>y> A,M,e), 


where the new variables A, y are A = i the functions f, g and h are as follows: 


= -3x1 - —y [y2) - . 

„ _ 76 3 ^ 

J2 = AxiX2 - -X2^ . 

/a = - 2 x 3 - ^2:35 . 

gi = -3eyi + VeX (x 2 )^ + V6g ( 2 / 2 )^ - ^^2/iS • 

76 3 „ 

92 = —^m2/2 - 2^2^ ■ 

hi = -TeA^Fi - l)xi ; Fi = 




(40) 

(41) 

(42) 

(43) 

(44) 

(45) 

(46) 

(47) 

(48) 


4.1. Slow-fast dynamical system 

The dynamical system dsni-gi]) is a singular perturbation system; specifically, it is a slow-fast system of first-order 
ordinary differential equations, where 2/1 is the fast variable [6ll4^ . Under the transformation N = eiV, the dynamical 
system (|591)-(|4T|1 becomes 

dx „ dA dyi dy2 dy 

dN~^^ ^ ’ 


( 49 ) 
















TABLE I: Fixed points analysis of the evolution equations ISSli-dlTj in the slow manifold for the exponential potential 


Point (xi,X 2 ,X 3 ) 

(yi,y2) 

P 

Eigenvalues 

Stability 


nr 

ncup 

We/.f 

'Wgup 

0 

(0,0,0) 

(0,0) 

h 

1 (-1,1,1,0) 

Unstable 

1 

0 

0 

0 

$ 

0* 

(0,0,0) 

( 2 / 1 , 0 ) 

0 

1 (-1,1,1,0) 

Unstable 

1 

0 

0 

0 

$ 

A± 

(±1,0,0) 

(0,0) 

p 


Unstable 

0 

0 

1 

1 

1 

Al 

(±1,0,0) 

(2/1,0) 

0 

(3,3,3t^,0) 

Unstable 

0 

0 

1 

1 

1 

R± 

(0,0, ±1) 

(0,0) 

p 

(-1,2, 2,0) 

Unstable 

0 

1 

0 

1 

3 

$ 

Rl 

(0,0, ±1) 

( 2 / 1 , 0 ) 

0 

(-1,2, 2,0) 

Unstable 

0 

1 

0 

1 

3 

$ 


which is a regular perturbation system The two dynamical systems (I39l)-6T|) and (|4^ are equivalent but in 
different “time scales”. 

In order to continue, we will use the Tikhonov’s theorem (^ . 1^ . and, we will study the fixed points of the 
dynamical system (|5ni) - (|^T|) in the slow manifold 52 (x, y, A,/i,0) = 0. This solution is called outer solution of the 
singular perturbed system. 

In the limit £ = 0, the dynamical system is equivalent to the algebraic dynamical system (EH), (gOl) and 
52 (x, y. A, 0) = 0; the last gives the slow manifold A {x2)^ + ^ ( 2 / 2 )^ = 0. Furthermore, when e = 0, from (I5S1) . we 
have that xx^ < 1; that is, for the variables x , x G [—1,1] holds. Moreover, since V {(/)) > 0, the parameter X2 is 
constrained as follows: X 2 G [0,1]. 

We consider the exponential potential V {4>) = Voe~'^^, where cr is a constant different from zero, then A = ct ; hence 
we have that ri = l, r 2 = ^. The fixed points of the outer solution with the physical value of the problem and the 
corresponding eigenvalues are given in table HI 

The fixed points of the table U exist for all values of the parameter A. Points O, O* correspond to the matter epoch, 
where the total EoS parameter is Wtot = 0. At the points A±, the universe is dominated by the kinetic energy 
of the scalar field = 1, V {<j)) = 0) , which means that the scalar field acts as a stiff fluid, i.e. = P(i)oa~^ , this 
solution corresponds to a non-accelerated universe. The fixed points i?±, correspond to the radiation epoch, where 
rij. = 1 and Wtot = Finally, the linearized system (EH, EH at the fixed points admits positive eigenvalues; that is, 
all the fixed points are unstable. 

However, in contrary to the “classical” minimal coupled scalar field cosmology, there are not the scalar field 
dominated and the tracker solutions [1^. We would like to note that, this analysis holds for very small values 
of the parameter AN = N — Nq; where Nq is the value which we consider as a starting point. 

When £ = 0, the regular perturbation system (H^ gives that the parameters x,X,y 2 :P are constants (that holds 
for any potential), and, the only non-constant equation is 


= V^A ( 3 : 2 )^ + ( 2 / 2 )^ = a = const, 

dN 


(50) 


that is, 2/1 (IV) = aN + 13, where a, /3 G K.. However, since £ = 0, the cosmological parameters are independent on 
2 / 1 ; that is, for very large values of the parameter {N — Nq) = In a, the cosmological parameters Hr, etc. will be 
constants. 

In order to understand better the effects of the GUP in the dark energy models, we perform numerical simulations 
for the dynamical system (EH'EH the power law model V {(j)) = Vo(j)~^ [1^, and, the hyperbolic model V (</>) = 
Vo (cosh (cr0) — I)*’ E^. 

In fig. [TJ we give the numerical solutions of the energy density evolution, and of the EoS parameter wgup for the 
the potential V {(f) = Vo(t>~^■ For the non-modified scalar field, the same plots can be found in [ 2 ^. We observe that 
in contrary with the “classical” quintessence scalar field, where in the tracking solution, the EoS parameter has the 
value W(j, = , whereas the wgup can cross the phantom divide line; however, the energy density evolution of the 

two models has similar behavior. 

Furthermore, in fig. [H we compare the evolution of the EoS parameters and for the modified EoS parameter 
WGUP for the hyperbolic model V {(j)) = Vo (cosh((T^) — 1)^. We observe that the scalar field modified by GUP can 
cross the phantom divide line, i.e. wgup < — 1- 

Finally, in order to understand the behavior of the slow-fast system (EH'EH: in hg- El we give the evolution of the 


parameter 


,rfyi 
' dN 


for the hyperbolic potential. From the fig. El we observe that for large AN parameter, 2/1 is fast, i.e. 

I can be measurable. For the numerical solutions, in figs, 
m El and El we selected £ = 10 which corresponds to the natural value of the parameter /3o , i-e. /3o — 1, However, 


2/1 ~ £ which means that the term eyi in equation 
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Energy density vs redshift in GUP for ^ 


0 1 2 3 4 5 6 7 

iog10(1+z) 


Evoiution of for V((t))-(^ ' 




FIG. 1: Numerical solution of the field equation in the dimensionless variables for minimally coupled scalar field cosmology 
modified by the GUP for the power law potential V (0) = The left fig. is the evolution of the energy density of the 

flm (solid line), fir (dashed line) and fl^ (dot line). For the numerical simulation we start from points in the radiation epoch, 
log 10 (1 + a) = 7.21. The right fig. is the evolution of the scalar field equation of state parameter (EoS) for the same model. 
For the numerical integration, we have considered /Iq = 1. 


we would like to remark that the behavior of the solutions is similar and for smaller values of the parameter e, i.e. 
/3o; since when the parameter e decrease, the fast variable yi becomes “faster”. 


5. CONCLUSIONS 


In this work, we have considered a modified quintessence scalar field Lagrangian which follows from the GUP. 
The modified Klein-Gordon equation for the scalar field admits a fourth-order perturbation term which follows from 
the quantum corrections. Furthermore, we have considered the gravitational action, and, with the use of Lagrange 
multiplier, we rewrote the gravitational action as that of a two scalar fields model with a mixed kinetic term; we remark 
that, the second field was introduced by the quantum corrections of GUP and it is not a real field. The proposed 
model, i.e. the action integral (1121) . it can be compared with that of a quintom model in which two fields take place. 
However, in our consideration there exist only one field which satisfies a fourth-order Klein-Gordon equation and in 
general it is different from the quintom models, moreover the coordinate transformation which diagonalize the kinetic 
term, which is defined by the field (p, and the Lagrange multiplier ip, is a singular transformation; which means that 
the inverse transformation can not be defined. 

Consider now the action integral of a quintom model with mixed kinetic term, which has been introduced in [s^ . 


Srn.ix, — 


d‘^3 






+ - f + Vd (</>) + U2 (V^) 


(51) 


where 0Oi V’o.o; are constants^, then in comparison with the action integral (I12L it can be seen that, ipQ = 0, V 2 (ip) is 
fixed, which follows from the Lagrange multiplier, and the constant, a, is related with the deformation parameter of 
the GUP, i.e. a ~ /3. Hence a quintom scenario with mixed terms in which the parameter, a, is small, can arise form 
the GUP. 

In the case of a spatially flat FLRW spacetime, we derived the field equations and we wrote the modified cosmological 
parameters. We show that the EoS parameter wgup for the scalar field has not a minimum boundary; that is, it 


In [Hi , 00 = bo = 1 
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Evolution of for V{(t))=V0{cosh{o$)-1 


Evolution of for V{(l))=V0(cosh(a(l))-1 




FIG. 2: Numerical solution of the EoS parameter and wgup for minimally coupled scalar field cosmology for the hyperbolic 
potential V {(f)) = Vq (cosh {(J(j>) — 1)^ with p = 0.2, and, a — 20 [^. The left fig. describes the evolution of the EoS parameter 
for the “classical” scalar field, whereas the right hg. is for the modified EoS parameter wgup- For the numerical simulation, 
we start from points in the radiation epoch, log 10 (1 + z) = 7.21, and, we consider j3o = 1. 


Evolution of the parameter |e dy,/dN| for the hyperbolic model 



FIG. 3: Numerical solution of the parameter edyi/dN (equation (1411) ) for the hyperbolic model of hg[2l We observe that, for 
large values of AN, the fast variable yi take values yi ~ and, when edyi/dN > 1, the wgup crosses the phantom divide 
line, which means that the term eyi of equation (13811 is measurable. For the numerical integration, we have considered /3o = 1. 
However, a similar behavior holds, and for different values of Pq. 


is possible to cross the phantom divide line, such as, wgup < ~1- Moreover, in order to study the cosmological 
evolution of the model, we wrote the field equations in the dimensionless variables, and, in the lapse time N = Ina. 
Since the dynamical system is a singular perturbation system we studied the existence of fixed points in the slow 
manifold. We have shown that, for an exponential potential, the outer solution of the field equations admits only 
unstable fixed points, which correspond to the matter dominated era, the radiation epoch, and, to the solution where 
the universe is dominated by the kinetic part of the scalar field. For “large” values of the lapse time, we demonstrate 
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that the cosmological parameters are constants. 

Finally, we have performed numerical simulations of the field equations for two well known potentials V (</>); in 
both cases, the EoS parameter wgup can take values lower than minus one. That result is important, since we can 
have a phantom behavior which follows for the quantum corrections of a quintessence scalar field model. 
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